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Introduction.  We  have  been  working  on  the  following  issues; 

1.  Numerical  methods  and  computational  experiments  for  the  multi-dimensional 

Inverse  Scattering  Problem  (ISP)  for  the  Helmholtz  equation  [3]. 

2.  Numerical  methods  for  the  3-dimensional  ISP  in  random  media  [6]. 

3.  Global  convergence  of  numerical  methods  for  3-D  ISP  with  time  dependent  data  [7]. 

4.  Phaaeless  1-dimensional  ISP  [8]. 

5.  Quasi-Reversibility  method  for  nonclassical  Cauchy  problems  [9]-[ll]. 

Overall  almost  all  our  approaches  axe  entirely  new.  In  peirticular,  we  note  that  our 
numerical  methods,  for  multidimensional  ISPs  see  [3],  [6],  [7]  and  [11]  work  directly  with 
PDEs  rather  than  with  corresponding  integral  equations.  And  this  is  a  crucial  aspect  in 
which  our  methods  differ  from  most  others.  This  report  is  supplemented  with  our 
reprints/preprints  [3],  [5],  [6]  [8]-[12].  For  this  reason,  we  discuss  all  the  issues  rather 
briefly. 

At  the  end  of  the  report  we  briefly  discuss  applications  of  our  results  to  ONR  ARI 
SEA  ICE  INITIATIVE  [2].  We  strongly  believe  that  our  research  is  very  well  suited  for 
this  initiative. 

1.  Numerical  Methods  and  Computational  Elxperiments  for  Multidimensional  ISP 

for  the  Helmholtz  Elquation  [3]. 

In  [13]  we  have  described  a  version  of  quasi-Newton  method  for  this  ISP  and 
presented  results  of  numerical  experiments.  We  have  noted,  however,  that  the  presence 
of  the  highly  oscillatory  function  in  Helmholtz  equation  negatively  affects  our 

computational  results  in  terms  of  both  computer  time  and  memory.  Here  the  unit 
vector  1/  represents  the  direction  of  the  incident  plane  wave. 

Recently  we  have  modified  our  method.  Namely,  let  u(x,  i/)  be  a  solution  of  the 
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forward  scattering  problem.  We  introduce  the  function  w(x,  «/)  =  u). 

The  reader  can  compare  Re[u(x,  i/)]  and  Re[w(x,  v)]  on  figures  1  and  2,  in  the  case  of 
the  scattering  from  a  circular  cylinder.  For  an  additional  discussion  see  [3].  Now  we 
can  use  a  smaller  number  of  Fourier  coefficients  of  the  function  w  to  obtain  the  same 
precision  as  before.  As  a  result,  our  modified  computational  procedure  requires  10  times 
less  com^Juter  time  than  before,  in  2-D  case.  For  instance,  now  the  reconstruction 
process  requires  8  sec.  of  CPU  time  instead  of  84  sec.,  see  [3]. 

Based  on  these  results,  we  hope  to  conduct  computational  experiments  in  the  3- 
Dimensional  case  in  the  near  future.  Clearly  such  results  are  important  for  ARI  [2]. 


2.  S-Dimensional  ISP  in  Random  Media  [6] 

We  have  noted  in  [6]  that  Forward  Scattering  Problems  (ISP)  in  the  random  media 
can  be  successfully  described  by  nonstationary  Schrbdinger  equation  (also  see  [17]).  To 
our  best  knowledge,  this  equation  has  never  been  used  by  mathematicians  to  study 
multidimensional  Inverse  Scattering  Problems  either  in  random  or  in  deterministic 
cases.  Physicists  accept  a  certain  “deterministic”  model  for  this  random  FSP,  see  [6], 
[17].  They  formally  derive  this  model  from  the  random  case.  In  this  model,  ensemble 
average  and  dispersion  of  the  random  density  are  engaged  as  real  and  imaginary  parts  of 
the  potential  accordingly.  We  have  developed  a  quasi-Newton  method  for 
“deterministic”  ISPs,  which  is  naturally  linked  with  the  above  mentioned 
“deterministic”  FSP.  We  have  established  the  convergence  rate  of  our  method.  This 
ISP  has  significant  applications  in  biology,  underwater  acoustic,  and  other  areas. 
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3.  Global  Convergence  of  Numerical  Method  for  3-Dimensional  ISP 
with  Time  Dependent  Data  [7] 

We  consider  the  following  potential  ISP  for  a  hyperbolic  equation.  Let  T  be 
positive  numbers  and  e  <  T.  Suppose  that  function  a(x)  e  C^(R^).  Let  this  function  be 
known  for  xe  {lxj<€}.  Determine  the  funtion  a(x)  for  xe  {/x/<^T}  assuming  that 

Ui^^=  Au+ a(x)u,  (x,  t)eR^x(0,  T), 

lxl=e  =  ‘P(^> 

where  the  function  p(x,  t)  is  known. 

In  [11]  we  have  developed  quasi-Newton  method  for  similar  ISP  assuming  that 
potential  a(x)  is  sufficiently  small.  Recently  we  obtained  a  fundamentally  new  result. 
Namely,  we  introduce  a  special  weighted  minimized  cost  functional  linked  with 
Carleman  estimates  see  [1],  [12].  Then  using  these  estimates  we  prove  that  our 
functional  is  strictly  convex  on  each  convex  compact  set  of  each  finite  dimensional  space 
(by  certain  choice  of  parameters  of  Carleman  estimate).  Thus  this  result  guarantees 
convergence  of  any  version  of  the  gradient  method,  if  the  starting  point  is  chosen 
appropriately.  W.  Symes  has  obtained  a  similar  result  for  the  one-dimensional  problem 
in  [18].  He  has  also  raised  the  question  about  the  3-D  case. 

4.  Phaseless  1-Dimensional  ISP  [8] 

In  the  classical  formulation  of  1-D  ISP  for  Schrodinger  equation,  one  has  to 
determine  the  potential  from  the  reflection  coefficient  R-(k).  In  many  important 
applications,  however,  only  the  amplitude  |R-(k)|  can  be  measured.  This  ISP  is 
tightly  connected  with  the  well  known  phase  problem  in  optics  see  [4],  [14],  and  [15]. 
For  the  fiirst  time,  we  have  proved  a  uniqueness  theorem  and  developed  rigorous 
numerical  method  for  this  phaseless  ISP. 
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5.  Quasi-Reversibility  Method  and  Stability  for  Nondassical  Cauchy  Problems  [5],  [9]- [11]. 

It  is  commonly  known  that  numerical  methods  for  the  nondassical  Cauchy 
problems  axe  important  value  for  numerical  methods  for  inverse  problems.  The  Quasi- 
Reversibility  (QR)  method  was  developed  in  [16].  But  convergence  rates  were  not 
established,  and  numerical  experiments  were  conducted  without  noise  in  the  data. 
Besides,  the  QR  method  has  not  been  derived  for  the  hyperbolic  timelike  Cauchy 
problem.  The  Caxleman  estimates  play  a  fundamental  role  in  the  QR  theory.  We  have 
shown  this  in  [5],  [9]-[ll].  Using  the  Caxleman  estimates  we  established  convergence 
rates  for  different  versions  of  the  QR  method.  Numerical  experiments  were  conducted 
both  in  the  case  of  the  Laplace  equation  and  the  wave  equation  with  noise  in  the  data. 
These  experiments  demonstrate  high  efficiency  and  stability  of  the  QR  method. 

6.  Applications  to  ONR  Sea  Ice  Initiative 

Caxefully  reading  [2]  we  have  concluded  that  our  research  is  closely  related  to 
challenges  of  the  ARI.  In  order  to  demonstrate  that,  we  draw  some  excerpts  from  [2] 
below  and  show  how  our  results  could  work  for  these  chedlenges. 

1.  [2,  p.  2]  “Sea  ice  is  inhomogeneous  and  highy  variable  at  all  spatial  scales  from 

millimeters  to  kilometers  ...” 

We  think  it  would  be  reasonable  to  assume  that  small-scale  variations  of  sea  ice  are 
random  and  large-scale  variations  aire  rather  deterministic.  Then  if  one  would  try  to 
determine  ensemble  averages  and  dispersions  of  corresponding  random  fields,  then  our 
results  [6]  will  work.  Besides,  we  hope  to  extend  results  [6]  on  Maxwell  systems. 
Further,  numerical  method  [6]  is  based  on  the  same  idea  as  ones  in  [3].  Therefore 
computational  experiments  conducted  in  [3]  could  help  to  make  computations  for  ISPs 
in  random  media.  Besides,  we  hope  to  conduct  numerical  experiments  in  the  3-D  case 
in  the  near  future,  and  this  CMe  is  definitely  more  realistic  than  2-D  cases. 
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2.  [2,  p.3]  “Can  correlation  functions  ...  be  successfully  computed  for  modelling 
purposes  ...” 

The  result  of  [6]  is  possibly  related  to  this  issue.  In  particular,  this  result  can  be 
applied  to  biological  media,  see  also  section  2  above.  All  our  results  concerning  3-D 
ISPs  are  closely  related  to  several  aspects  metioned  in  the  other  questions  in  Section  D 
of  the  ARI. 

3.  [2,  p.3],  bottom  and  top  of  p.  4 

“There  is  need  for  the  development  and  validation  of  time- dependent  models  for 
beam  spread  in  sea  ice  ..  .  Mathematical  models  for  EM  inverse 
scattering  must  be  developed  to  address  these  needs.  ” 

Our  results  in  [7]  and  [11]  respond  positively  to  these  challenges.  The  method  in  [7] 
is  of  special  value  since  it  is  a  global  method,  i.e.,  a  good  first  guess  is  not  required. 
Furthermore  we  hope  to  extend  our  methods  to  the  more  difficult  and  useful  equation 
'itt  =  a^(x)Au. 

4.  [2,  p.7]  “.  .  .  have  little  theoretical  guidance  as  to  which  measurements  are  most 
essential  and  the  measurement  accuracy  required  to  draw  strong  conclusions” 

All  our  results  together  with  uniqueness  eind  stability  theorems  [l],  [9]  -  [12]  can 
provide  a  strong  theoretical  basis  for  this  challenge.  Besides,  since  numerical 
experiments  show  high  stability  of  our  methods,  see  [3],  [9]  Jind  [10],  they  could  help  to 
predict  the  accuracy  required.  Furthermore,  there  are  a  large  number  of  measurements 
for  which  only  the  amplitude  of  the  signal  is  available,  but  not  the  phase  [4].  Thus  the 
phase  problem  in  optics  arises  very  naturally  in  these  situations.  Therefore  our  results 
[8],  [14],  and  [15]  can  cert£iinly  be  of  use  in  such  cases. 
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Figure  Captions 

1.  Re[u(x,  j/)],  where  v  =  {\,  0),  and  u  is  scattering  field  from  a  circular  cylinder  with 
the  radius  a  =  e  =  0.03  inside  of  the  cylinder,  t  =  0  outside  of  the  cylinder  and 

k  =  5.3.  Here  e  is  a  conductivity,  see  [3]. 

2.  Re[u(x,  with  the  same  coefficient  e  as  on  the  Figure  1.  The  features 

of  the  field  are  much  more  clear  than  ones  on  Figure  1. 
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